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We present a next-to-leading order (NLO) QCD analysis of 
unpolarized and polarized deeply virtual Compton scattering 
(DVCS) amplitudes, for two different input scenarios, in the 
MS scheme. We illustrate and discuss the size of the NLO 
effects and the behavior of the amplitudes in skewedness, £, 
and photon virtuality, Q 2 . In the unpolarized case, at fixed 
Q 2 , we find a remarkable effective power-law behaviour in f, 
akin to Regge factorization, over several orders of magnitude 
in £. We also quantify the ratio of real to imaginary parts of 
the DVCS amplitudes and their sensitivity to changes of the 
factorization scale. 



I. INTRODUCTION 

Deeply virtual Compton scattering (DVCS) ]!]-[)), 
7* (?) + P{P) —* lil') + P{P')i i s the most promis- 
ing process for accessing generalized parton distributions 
(GPDs) which carry new information about 

the dynamical degrees of freedom inside a nucleon. GPDs 
are an extension of the well-known parton distribution 
functions (PDFs) appearing in inclusive processes and 
arc defined as the Fourier transform of non-local light- 
cone operators sandwiched between nucleon states of dif- 
ferent momenta, commensurate with a finite momentum 
transfer in the t-channcl. These distributions are true 
two-particle correlation functions and contain, in addi- 
tion to the usual PDF-type information residing in the 
DGLAP jl3) region, supplementary information about 
the distribution amplitudes of virtual meson-like states 
in the nucleon in the ERBL jTJJ] region. 

We recently presented a full numerical solution of the 
associated renormalization group equations at next-to- 
leading order (NLO) accuracy, for unpolarized and po- 
larized distributions, using realistic input models jL5|, as 
well as a complete NLO QCD analysis of DVCS observ- 
ables Jlq| . To achieve this we calculated the real and 
imaginary parts of unpolarized and polarized DVCS am- 
plitudes, T^vqs, at NLO which are related to the triple 
differential cross section on the lepton level (which also 
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includes the Bethe-Heitler (BH) process) via 
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and M is the proton mass. The dependent variables 
Q 2 = —q 2 ,Xf,j = Q 2 /2P ■ q,t and <fi are minus the pho- 
ton virtuality, Bjorken x, the momentum transfer to the 
proton squared and the relative angle between the lepton 
and proton scattering planes Jl7t , respectively. Variable 
y = Q 2 /xbjS, where S is the total lepton-proton center 
of mass energy. For the DVCS process in the proton rest 
frame, y is the fraction of the energy of the incoming 
lepton carried by the photon. 

The real and imaginary parts of the polarized and un- 
polarized DVCS amplitudes are interesting in their own 
right, since each can be independently accessed experi- 
mentally via various asymmetries which exploit the <fi- 
dependence of the interference term Thus the 

predictions that we give here for the real and imaginary 
parts of the amplitudes at NLO, their behavior in the 
skewedness, £ = Xbj, and Q 2 , as well as the relationship 
between real and imaginary parts, which all depend on 
the GPDs, can be directly tested by experiment. A NLO 
analysis of the DVCS amplitudes for large C was carried 
out in fl9|| and, at the limited points where a comparison 
is possible, we agree with their results. 

This paper is structured as follows. In Section [n] we 
define the input model for the GPDs and briefly review 
their NLO evolution. We state the necessary factorized 
convolution integrals and explain their exact technical 
implementation, via subtractions. In Section [II we give 
detailed NLO results for the real and imaginary parts 
of the unpolarized (subsection III A ) and polarized (sub- 
section IIIB) DVCS amplitudes, comparing them to LO 



results using the same input GPDs and commenting on 
their sensitivity to the input GPDs and the factorization 
scale, fj, 2 . We also give the Q 2 and C dependence of the 
ratio of real to imaginary parts and quantise the impor- 
tance of the ERBL region to the real part. Wc close our 
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discussion in Section |y| with a statement on the gen- 
eral structure of NLO and NNLO corrections and briefly 
conclude in Section For convenience in an appendix 
(se ction |vT[ ) we restate the NLO coefficient functions |l|] 
(in VI AD and give analytic results for their integrals (in 



VI B), which are required to implement the subtractions 
in section 



HI. 



II. 



FACTORIZATION THEOREM AND 
DEFINITIONS 



There are many representations for GPDs in the litera- 
ture [12|]2C| l . We chose to work in a particular rep- 
resentation identical to the non-diagonal representation 
defined in po|] which is a natural one when comparing to 
experiments. We use 



T s ^' v / A {X,C^\t) = 
H a < v / A (v, C, ii 2 , t) t H a ^ A (- Vl C, m 2 , t) 



(l-C/2) 

H3- v / A (v,(,^,t)±H9> v / A (-v,t,^,t) 



(l-C/2) 
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where v = (X - C/2)/(l - C/2), S(a) and g refer to the 
quark singlet for flavour a and the gluon, respectively, 
V and A stand for unpolarized (vector) and polarized 
(axial- vector) cases, taking the upper and lower signs, re- 
spectively. This representation is different from the usual 
one (see for example [To|]) which is defined symmetrically 
with respect to the incoming and outgoing nucleon plus 
momentum (defined on the interval v G [—1,1] and sym- 
metric about v = 0). The GPDs in eq. (p|) have plus 
momentum fractions (on the interval X G [0, 1]) with re- 
spect to the incoming nucleon momentum, P, in analogy 
to the PDFs of inclusive reactions, with the ERBL re- 
gion in the interval X G [0, C] and the DGLAP region 
in the interval X G [Cjl]- The transformation between 
the symmetric and non-diagonal representation is given 
in p0| . Furthermore, within the non-diagonal represen- 
tation £ = xi,j = —q 2 /2P ■ q and the symmetry of the 
GPDs, which was previously manifest about v — 0, is 
now manifest about the point X = £/2. 

We build input distributions, F s{a ^ 9 (X, C, Ql), at the 
input scale, Qo, with the correct symmetries and proper- 
ties, from conventional PDFs in the DGLAP region, for 
both the unpolarized and polarized cases, by employing 
the factorized model due to Radyushkin |2]}| , which is 
based on double distributions. The H functions (sym- 
metric GPDs) required for eq.(||) are related to the lat- 
ter via the following reduction formula (factoring out the 



overall ^-dependence): 

H(v, C) = J 1 dx> J 1 '* ' dy'S f x ' + ^- v) F(x', y') . 

(3) 



The double distributions, F l ' V t A , are a product of a 
profile function, 7r l , and a conventional PDF, f l ' V ^ A , (i ~ 

F"( a \x',y')=^(x',y')f^(x') 
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where 



f 9 (x) = xg(x, Qo)Q{x) + \x\g(\x\,Q )Q{-x) 
f« a \x) - q a (x, Q )Q(x) - (q a )(\x\, Q )Q(-x) . (5) 

The profile functions are chosen to guarantee the correct 
symmetry properties in the ERBL region which are pre- 
served under evolution, as we explicitly illustrated in [jl5| . 
Their normalization is specified by demanding that the 
conventional distributions are reproduced in the forward 
limit at the input scale: F l {X, ( -> 0, Q ) f(X, Q ). 

In addition to the contributions from the double dis- 
tributions the unpolarized singlet GPDs also contain a 
so-called "D-term" [^2|,^3| , which is only non-zero in the 
ERBL region, and ensures the correct polynomiality in 
C ]23|] of the moments in X of the GPDs. There is an 
equivalent term in the unpolarized gluon distribution but 
apart from its symmetry nothing is known about this 
function, thus we chose to set it to zero. 

Within the above class of input model, we specify 
two particular input models for the GPDs by using two 
sets of inclusive unpolarized/polarizcd PDFs (for use in 



eqs.(|i|,|)), i-c. GRV98 / GRS V00 @] with A, 



246 MeV and MRSA'/GS(A) 



with A 



(4,NLO) 

QCD 

(4,NLO) 

QCD 



231 MeV at the common input scale Qq = 4 GeV 
174 MeV for both sets. Using two dif- 



and A 



(4,LO) 

QCD 



ferent choices allows us to investigate the sensitivity of 
the amplitudes to the choice of input. The GPDs are 
then evolved in LO and NLO using our newly developed 
evolution code |l5| . Note that there are two sets of evo- 
lution equations which have to be solved simultaneously, 
one for the ERBL region Q and one for the DGLAP 
]l3] | region, with the ERBL one being dependent on the 
evolution in the DGLAP region, whereas the evolution in 
the DGLAP region is independent from the evolution in 
the ERBL region (for more details on the NLO skewed 
evolution and NLO coefficient functions see |l5 26|). 
The factorization theorem pJW| proves that the DVCS 



amplitude takes the following factorized form (in the 
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non-diagonal representation) up to terms suppressed by 
0(1/0): 



2-C 



c 
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In the following we will initially set the factorization 
scale, fi 2 , equal to the photon virtuality, Q 2 , and will later 
investigate its variation. Henceforth, we will suppress 
the factorized i-dependence since all of our predictions 
will be made for t = 0. The LO and NLO coefficient 
functions, T l * y / A , are taken from eqs. (1 4-17) of [|l9| and 
are summarized in the appendix VIA . They can have 
both real and imaginary parts, depending on the region 
of integration, which in turn generates real and imaginary 
parts of the DVCS amplitudes. For the integrals over 
the range X e [0,1], the coefficient functions contain 
singularities associated with the point X = C, which are 
regulated using a "+ie" prescription. For our numerical 
integration, we choose to regulate these integrals using 
the Cauchy principal value prescription (denoted P.V.) 
as follows: 



2X 

T 7 



P.V. I >1X T[ — -1 )t(x,c,Q 2 ) 



^ dXT(^- lj (F(X, c, Q 2 ) - HC, C, Q 2 )) 



(7) 



Each term in eq. (^) is now either separately finite or 
only contains an integrable logarithmic singularity. This 
algorithm closely resembles the implementation of the + 
rcgularization in the evolution of PDFs and GPDs. Note 
that the first integral in eq. (fjj) (in the ERBL region) is 
strictly real, however, the second and third terms contain 
both real and imaginary parts (which are generated in the 
DGLAP region). This definition leads to the following 
formulas for the real and imaginary parts of the DVCS 
amplitudes: 



f l dX T S ^> V ' A - 1, F S W A {X, C, M 2 , t) T f° 

1 dX T S ^' V ' A J*MW(X, C, M 2 , t) 



ReT^(C0 2 ) = E e a C 2 -^) 
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(9) 



where z = 2AT/C — 1 and for convenience we have sup- 
pressed the scale and explicit quark flavour dependence 
of the GPDs on the right hand sides. The real and imag- 
inary parts of the unpolarized and polarized DVCS am- 
plitudes were computed using a FORTRAN code based 
on numerical integration routines (for more details see 
p7|). We implemented the exact solution to the RGE 
equation for a s , in LO and NLO as appropriate, in our 
calculation to be consistent throughout our analysis. 
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III. NLO AND LO DVCS AMPLITUDES 

In the next two subsections we present results for the 
real and imaginary parts of the NLO DVCS amplitudes 
for the unpolarized and polarized cases, comparing them 
to LO results (using the same input GPDs) throughout. 
This allows us to quantify the effect of moving from LO 
to NLO. We first plot the absolute values, then the ratio 
of real to imaginary parts (in Q 2 for fixed £ and in £ 
for fixed Q 2 )- Finally, we discuss the influence of the 
ERBL region on the real part of the amplitudes and the 
factorization scale dependence. 



from small to large £ reflecting the strong decrease in 
the GPDs. They generally increase with increasing Q 2 , 
albeit moderately, reflecting the expected ln(Q 2 ) behav- 
ior. In fact, an approximate scaling is observed at LO 
and NLO, but only sets in at large Q 2 for small £ in the 
imaginary part at NLO. 

In Fig.([l]) and subsequent figures a comparison of the 
LO and NLO curves at the input scale Q 2 — 4 GeV 2 
reveals the effect of including NLO coefficient functions 
alone. Of course at higher Q 2 NLO effects are included 
in both the coefficient functions and in the evolution of 
the GPDs. For the numerical size of the NLO corrections 
in the evolution of the GPDs alone see jL5|. 



A. The unpolarized DVCS amplitudes 



Unpol. Re and Im gluon DVCS Amplitude in NLO 



In Fig.© we plot the real and imaginary parts of 
the unpolarized quark singlet amplitude (cf. eq. (||)) at 
LO and NLO. We observe that MRS A' and GRV98 give 
broadly similar results. 



Unpol. Re/lm singlet DVCS Amplitude in LO/NLO 
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FIG. 1. The Q 2 -dependence of the real and imaginary 
parts of the quark singlet DVCS amplitude. The solid 
(dashed) curve is the real part in LO (NLO) and the dotted 
(dashed-dotted) curve is the imaginary part in LO (NLO). 

The NLO corrections are generally fairly large and tend 
to decrease the value of the real part and increase the 
value of the imaginary part of the quark singlet ampli- 
tude (bearing in mind that the NLO imaginary amplitude 
is lower at the input scale, due to the inclusion of the 
NLO coefficient function). The relative correction, R, in 
moving from LO to NLO (i.e. R = (NLO-LO)/LO)) can 
be inferred from Fig.©. R is typically -20% to -40% 
for the real part for both values of £. For the imaginary 
part at £ = 0.1 the corrections are moderate (\R\ < 20%), 
whereas at small £, R can be as large as +50% for GRV98 
at large Q 2 . The amplitudes drop dramatically in going 
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FIG. 2. The Q 2 -dependence of the real (solid line) and 
imaginary (dotted line) parts of the unpolarized gluon DVCS 
amplitude, for two representative values of £. 

In Fig. (^) we show the real and imaginary parts of the 
unpolarized gluon amplitude (cf. eq. (^J)) which starts 
at NLO. We note firstly that the gluon contribution is 
of the same order of magnitude as the quark singlet one, 
although it is suppressed by a s /2n and secondly that 
both the real and imaginary parts are large and nega- 
tive. This explains the strong variation of the azimuthal 
angle asymmetry (AAA), in moving from LO to NLO, ob- 
served in |l6| (for the case of a dominant Bethe-Heitler 
contribution) which is directly proportional to the real 
part of a combination of DVCS amplitudes. Otherwise, 
the gluon mirrors the NLO quark singlet amplitudes in 
its behavior in Q 2 and for large and small £. 

Of course for the physical amplitude at NLO one must 
add the quark singlet and gluon contributions together. 
For the real part at small £ = 10~ 4 this leads to a rela- 
tive relative change in moving from LO to NLO of 40% 
(MRSA) and 80% (GRV98) at the input scale to 75% 
(MRSA) and 80% (GRV98) at Q 2 = 100 GeV 2 . For 
large £ = 0.1 the relative change in the real part at both 
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scales is about 60%. For the imaginary part at small 
C = 10~ 4 the relative changes for MRSA' are 35% (at the 
input) and 15% (at Q 2 = 100 GeV 2 ). For GRV98 one 
finds 50% (input scale) and 2% (evolved scale) changes. 
For the imaginary part at large £ = 0.1, for both MRSA' 
and GRV98, one finds about 30% (input scale) and 10% 
(evolved scale) changes. Note that the NLO changes in 
the physical amplitude decrease as Q 2 increases in line 
with the expectation of perturbative QCD that the NLO 
corrections should die out as Q 2 — > oo. 

In Fig. (|3|) we show the ratio of the real to imaginary 
parts for both the quark singlet and the gluon amplitudes 
as a function of Q 2 . We note that for small £ the ratios 
can be as large as 45% (a similar value was found for the 
closely related process of high energy J/ ip photoproduc- 
tion in |^8| ) . The quark ratios are slightly different for the 
two inputs. The greatest contrast is seen at small £: the 
NLO case is basically flat in Q 2 which differs markedly 
from LO which rises with Q 2 . Both cases are similar at 
large The gluon ratio is remarkably similar for the two 
input sets, at both small and large £, for the whole Q 2 
range considered. 



Ratio of unpol. Re/lm DVCS Amplitude in LO/NLO 
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FIG. 3. The ratio of real to imaginary parts of the unpolar- 
ized quark singlet and gluon DVCS amplitudes, as a function 
of Q 2 . The solid (dotted) curve is for the quark singlet in LO 
(NLO) and the dashed curve is for the gluon in NLO. 

We now turn to the ^-dependence for fixed Q 2 which 
is shown in fig.(|4|) for the quark singlet and fig.(^) for 
the gluon. The most striking feature for the quark 
singlet case is that the amplitudes exhibit an effective 
powerdike behavior in £ over basically the whole range 
(C £ [0.0001,0.2]), as illustrated by the straight lines in 
fig.(Q). A simple two parameter fit of the type aoC A ° 
works remarkably well up to about £ = 0.1. The best 
fit is obtained with a four parameter fit, of the type 
ci£ Al (1 + C2C A2 ), which can reproduce most of the curves 



on the few percent level, with Ao and Ai within 5% — 20% 
of each other, and A2 small. The simple two parameter 
fit works best for the imaginary part of the amplitudes 
where we obtain a value for Ao between —1.1 and —1.25 
with a moderate growth in Q 2 as expected from measure- 
ments of the slope of F2 and HERA diffractive processes. 
For the real part of the amplitude the simple fit starts to 
decrease in quality around £ = 0.05 (depending on the 
input). 

Similar power-like behaviors have been observed in 
many small- a; processes at at HERA, including diffrac- 
tive DIS in which the x ^-dependence of diffractive struc- 
ture functions is known to factorize for small xp. This 
behavior, known as Regge factorization, is predicted for 
high energy processes within Regge theory given the pos- 
tulate of a universally exchanged Regge pole, known as 
the Pomeron. Since the relationship between the phe- 
nomenological Regge Theory and perturbative QCD re- 
mains unresolved, we find the observation of this single 
power in our numerical perturbative QCD calculation re- 
markable and very interesting. 



Unpol. Re and Im singlet DVCS Amplitude in LO/NLO 
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FIG. 4. The real and imaginary parts of the unpolarized 
quark singlet DVCS amplitude, as a function of £. The solid 
(dashed) curve is the real part in LO (NLO) and the dotted 
(dashed-dotted) curve is the imaginary part in LO (NLO). A 
remarkably simple behaviour is observed in £ which is close 
to a single power over a wide range. Because Q 2 is fixed, 
this behaviour in £ translates directly into a single power in 
energy. 

Note that we do not claim to have derived this power- 
like behaviour from first principles. The a nalyti c forms 
for the coefficients functions (see appendix VIA ) would 
seem to favor a more complicated sum of logs in £, par- 
ticularly when the interplay with the effective behaviour 
in £ of the GPDs, at a given Q 2 , is taken into account. 
Hence, the fact that a single power apparently works for 
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such a large range in £ for DVCS is somewhat surprising. 
Naively one may expect that a sum of two or more powers 
would be required to give a reasonable fit. Perhaps this 
indicates that DVCS always proceeds through partonic 
configurations which lie in the same universality class, i.e. 
are self-similar. The reason why these self-similar config- 
uration seem to be of importance beyond the "diffractive 
region" can be understood if one examines the behavior 
of the integrand in the convolution integrals in eqs. (||, 
|9|) . One observes that a large contribution to the integral 
and ultimately to the imaginary part of the amplitude 
itself stems from the region in X very close to £, even 
for larger values of £, leading to the self-similar behavior. 
This is not too surprising given the steep rise of the GPD 
in the DGLAP region towards C and that the coefficient 
functions are singular at £, where even the implementa- 
tion of the principal value in Sec. ( |jl| ) leaves an integrable 
singularity. 

This situation is somewhat altered for the real part 
of the amplitude where the ERBL region plays a very 
important part, as we will see. There, although the coef- 
ficient function is singular, the symmetry of the GPD in 
the ERBL region makes the value of the integral some- 
what less dependent on the region near £, i.e. less singu- 
lar, especially for large £. Although the value of the real 
part also depends on an integral over the DGLAP region 
which contributes more to the power like behavior, this 
dependence progressively decreases as Q grows. 

The physical picture which is emerging is the follow- 
ing: the region near £ corresponds to large light-cone dis- 
tances for the operators which translates into two parton 
configurations, either a qq pair in the ERBL region or 
a quark leaving and then returning to the proton in the 
DGLAP region, where one of the partons, either the q 
or the returning q, carries virtually zero momentum frac- 
tion, X — C, in either the + or — direction on the light 
cone and the other carries approximately £, i.e. very 
asymmetric configurations seem to have a disproportion- 
ate weighting in the amplitude. The coefficient functions, 
with their singularity structure at £, weight these config- 
urations much more heavily in the amplitude than other, 
more symmetric, configurations. 

Although asymmetric configurations become rare as 
one approaches the valence region of large £, because 
they should be mainly found in the sea, they are still 
enhanced by the singularity structure of the coefficient 
functions. This then leads us to the conclusion that, 
although the sea is small at large £, DVCS still proceeds 
largely through sea configurations in the valence region, 
thus its relative scarcity compared to DIS in the valence 
region. Physically, at large X w £, one has to strike an 
unusual fluctuation in the proton to emit a real photon, 
while still leaving the proton intact. 



Unpol. Re and Im gluon DVCS Amplitude in NLO 
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FIG. 5. The modulus of the real and imaginary parts of 
the unpolarized gluon DVCS amplitude as functions of for 
fixed Q 2 . The solid curve is the modulus of the real part and 
the dotted curve is the modulus of the imaginary part of the 
gluon amplitude. 

For the gluon contribution, shown in Fig.(^[), we find 
the same behavior as in the case of the quark singlet. 
Note that we took the modulus of both the real and 
imaginary parts, since they are actually negative, in or- 
der to produce a log- log plot. Performing the same type 
of fits as in the quark case, one obtains similar num- 
bers for Ao between —1.14 and —1.28 in the two param- 
eter fit and again very similar ones in the four parame- 
ter fit (5% — 25% variation in the powers). The quality 
of the two parameter fit starts to decrease rapidly for 
a ( ~ 0.05. Nevertheless, the explanation given in the 
quark singlet case is still applicable in the case of the 
gluon. 

In Fig. (||) we show the ratio of real to imaginary parts 
in £ for fixed Q 2 . Again we note the remarkable similarity 
between the gluon curves for both inputs, both in shape 
as well as absolute values. A very mild growth is seen 
for the quark singlet ratio in £, whereas the gluon varies 
quite strongly in £. 

Using Fig. (||) we can make a simple test to check 
whether we have computed the real part of the amplitude 
at small £ properly by employing a dispersion relation for 
the unpolarized amplitudes at small £ (for more details 
see p8[): 



Rc Tpycs 
Im Tdvcs 



— tan 



- 1 



(10) 



For the fitted values of A we obtain a ratio which is in very 
good agreement with the values in Fig. (^) up to about 
( = 5x 10~ 3 for both the quark singlet and the gluon. 
This confirms the self-consistency of our calculation. 
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FIG. 6. The ratio of real to imaginary parts of the unpolar- 
ized quark singlet and gluon DVCS amplitude, as a function 
of f, at fixed Q 2 . The solid (dotted) curve is the ratio in 
LO (NLO) and the dashed curve is the ratio for the gluon in 
NLO. 

Next we discuss the relative importance of the ERBL 
region to the value of the real part of the DVCS am- 
plitude, starting with the quark singlet. On inspec- 
tion of the relative contribution of the ERBL integrals 

(X G [0,C]) in eqs. (§, §) we nnd that at sma11 C the 
ERBL region integral has a relative contribution between 
90% at the input scale and 140% at Q 2 = 100 GeV 2 
(100% and 50% respectively in LO) of the value of the 
amplitude, i.e. there is a large cancellation between the 
subtraction term and the X € [£, 1] integral with both 
of them being substantially larger, individually than the 
[0, C] integral. As one increases £ the relative importance 
of the ERBL region drops to 50% at the input scale and 
130% at Q 2 = 100 GeV 2 (80% and 30% respectively in 
LO), however, the subtraction term now starts to dom- 
inate the value of the amplitude. This observation is in 
line with our previous argument of the importance of very 
asymmetric parton configurations. Remember, that the 
subtraction term is directly proportional to the GPD at 
(. Also note that going from LO to NLO seems to change 
the relative importance of the ERBL region, especially its 
Q 2 behavior. 

Turning now to the gluon we make a slightly different 
observation. Firstly, at small ( the relative contribution 
varies from 40% — 60% in going from the input scale to 
Q 2 = 100 GeV 2 again there is a large cancellation be- 
tween the previously mentioned terms, but the change 
in Q 2 is not very dramatic. Increasing £ one finds an 
increase of the relative importance to 60% — 80% in go- 
ing from the input scale to our large Q 2 value. Again 
the subtraction term becomes more and more dominant, 
relatively speaking, and thus our interpretation from the 



quark case carries over to the gluon case. 

To close our analysis of the unpolarized case, we discuss 
the scale dependence of the DVCS amplitudes. We varied 
the factorization scale, /z 2 , from Q 2 , used above, to Q 2 /2 
and 2Q 2 for both sets and found the following variations, 
where the two sets agree fairly well with one another. At 
small C, we found a small variation at the input scale of 
3 — 5% which increases to about 15% — 30% at large Q 2 
(we used Q 2 = 100 GeV 2 , for the large scale, throughout) 
for both real and imaginary part of the quark singlet. For 
the gluon we find a variation of about 15% for the real 
and imaginary part at the input scale which reduces to 
about 10% for the real part and to about 2% for the 
imaginary part. At large £, we find similar variations 
as the authors of jy|, i.e. around 5% for both the real 
and imaginary parts of the quark singlet and around 10% 
for the gluon, for both the input scale and at large Q 2 . 
In summary one can say that the scale dependence is 
not troublesome and that the uncertainties due to the 
chosen GPD are still much larger than those due to the 
factorization scale. 



B. The polarized case 

For the polarized case we proceed exactly as we did for 
the unpolarized case starting in Fig. (R) with the absolute 
values of the real and imaginary parts of the polarized 
quark singlet amplitudes as functions of Q 2 . Again the 
relative impact of the NLO corrections to the quark sin- 
glet can be inferred from this figure. One immediately 
notices for small £ the very large NLO corrections for 
the GS(A) input (\R\ can be as large as 70%), whereas 
for the GRSV00 input ('standard scenario', i.e. unbro- 
ken sea) the corrections are more moderate (\R\ < 20%). 
The very large corrections can be easily explained since in 
|]l5| it was shown that the GPD evolution drastically al- 
ters the shape of the GS(A) distribution in fact inverting 
its shape, whereas the shape of GRSV00 was almost un- 
changed and the absolute value changed only moderately. 
The difference in shape at small £ was due to radically 
different assumptions about the polarized sea distribu- 
tion. At large £, where the shape of the two input sets is 
similar we find that the shape of the amplitudes in Q 2 is 
also very similar, whereas their absolute values differ. 
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FIG. 7. Real and imaginary parts of the polarized quark 
singlet DVCS amplitude, as a function of Q 2 , for fixed £. 
The solid (dashed) curve is the real part in LO (NLO) and 
the dotted (dashed-dotted) curve is the imaginary part in LO 
(NLO). 

The real and imaginary parts of the polarised gluon 
amplitude are plotted in Fig. (||). Note that the real and 
imaginary parts of the amplitude are positive for small £, 
exactly the opposite to the unpolarized case. For large £, 
the real parts are positive, again in contrast to the unpo- 
larised case. The imaginary part at large C starts positive 
but becomes negative under evolution. This behavior is 
due to the particular shape of the polarized gluon GPD 
at small and large ( (sec fl5|). 
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FIG. 8. The real and imaginary parts of the polarized gluon 
DVCS amplitude, as a function of Q 2 , for fixed The solid 
curves show the real part and the dotted curves show the 
imaginary part of the gluon amplitude. 



In Fig. (|9|) we show the ratio of real to imaginary parts 
as a function of Q 2 , for fixed £• At large £, we again ob- 
serve flattening at large Q 2 : the GRSVOO result is flat- 
ter than the GS(A) case, especially at NLO, where we 
still observe strong variations in Q 2 . For large Q 2 , the 
gluon ratio is very similar for both sets and fairly flat in 
Q 2 . We did not plot the gluon ratio at large £ since it 
has a very large value near the point where the imagi- 
nary part changes sign and thus would have completely 
swamped the quark result which we find more interesting 
here. Note that our LO ratio for the GS(A) model is in 
agreement with the values obtained in [ p"8"| ]. 

In Fig. ( |Io|) we plot the ^-behaviour at fixed Q 2 and 
find a single power-like behaviour only for very small (, 
(up to about 3 x 10 -3 ), seen explicitly on the log-log plot 
for GS(A). Performing the same type of fits as in the un- 
polarized case, i.e. a simple two and four parameter fits, 
reveals a Ao between —0.4 to —0.55, and a similar story 
for the four parameter fits. However, the second power, 
02, is now substantially larger than in the unpolarized 
case, in order to be able to describe the large ( behavior. 
The four parameter fit is able to describe the behavior of 
the amplitudes on the few percent level. Thus the sim- 
ple single-power Regge-type effective behavior observed 
in the unpolarized case is only valid at small £ in the 
polarized case. This is because the polarized sea dies out 
even more quickly with increasing £, than the unpolar- 
ized one: therefore the highly asymmetric configurations 
necessary to produce a universal behavior in £ become 
very rare. At large £ they are almost completely gone. 
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FIG. 9. The ratio of real to imaginary parts of the polarized 
quark singlet and gluon DVCS amplitudes, as a function of 
Q 2 , at fixed 5- The solid (dotted) curves show the ratio for 
the quark singlet in LO (NLO) and the dashed curves show 
the ratio for the gluon in NLO. 
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FIG. 10. Polarized real and imaginary parts of the quark 
singlet DVCS amplitude, as a function of £, for fixed Q 2 . 
The solid (dashed) curve is the real part in LO (NLO) and 
the dotted (dashed-dotted) curve is the imaginary part in LO 
(NLO). 



Turning to the gluon in Fig. ( |11| ) we illustrate that the 
behavior in £ at an evolved scale, Q 2 = 9 GeV 2 , is very 
similar in shape and size for the two input sets, despite 
the fact that they start off very different at the input 
scale, Q 2 = 4 GeV 2 . The evolution forces the distribu- 
tions to be quite similar very quickly. 

Pol. Re and Im gluon DVCS Amplitude In NLO 




FIG. 11. Real and imaginary parts of the polarized gluon 
DVCS amplitude, as a function of £, at fixed Q 2 . The solid 
curve is the real part in NLO and the dotted curve is the 
imaginary part. 

For the GS(A) input, at small (, the relative changes 
in going from LO to NLO in the physical polarized am- 
plitude, i.e. the sum of polarized quark singlet and po- 



larized gluon, are found to be about 15% for both the 
real and imaginary parts at the input scale and about 
8% for the imaginary part and 40% for the real part at 
Q 2 = 100 GeV 2 . At large C = 0.1, the changes in the 
real part vary from about 90% at the input scale to 25% 
at Q 2 = 100 GeV 2 and those in the imaginary part vary 
from about 25% (at the input scale) to 15% at the evolved 
scale. For GRSVOO we also find a decrease in the varia- 
tion which at small £ drops from 60% at the input scale 
to 25% at Q 2 — 100 GeV 2 for the imaginary part, and 
from 35% to 5% for the real part. At large ( there is a 
decrease from about 25%, at the input for both real and 
imaginary parts, to about 2% and 15% for the real and 
imaginary parts, respectively, at the evolved scale. As in 
the unpolarized case the NLO corrections to the physical 
amplitude decrease as Q 2 increases except in the case of 
the real part at small ( for GS(A) which is in line with 
the dramatic NLO evolution effects in the ERBL region. 

The relative change induced by NLO corrections of the 
quark singlet can be inferred directly from Fig. [l^. 
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FIG. 12. The ratio of real to imaginary parts of the polar- 
ized quark singlet and gluon DVCS amplitudes, as a function 
of ( at fixed Q 2 . The solid (dotted) curve is the ratio in LO 
(NLO) and the dashed curve is the ratio for the gluon, which 
starts at NLO. 

The ^-dependence of the ratio of real to imaginary 
parts, plotted in Fig. (p"2|), is rather different to the unpo- 
larized case, with a strong fluctuation in the gluon due to 
a sign change of one of the amplitudes around the spike 
(we sample at a finite number of points in £). The real 
part is typically much larger then the imaginary part, in 
contrast to the unpolarized case. 

Concerning the issue of the importance of the ERBL 
region, we find a similar picture to the unpolarized case, 
although the relative percentage contribution is smaller 
by about a factor of two, for both quarks and gluons. 
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The influence of the subtraction term on the value of the 
real part of the amplitude is considerably less than in 
the unpolarized case, thus the ERBL region is weighted 
more heavily in the final value of the amplitude. This 
is in line with our observation of the deviations of the 
amplitudes from a single power in (, for relatively small 
values of £, compared to the unpolarized case (for which 
the ERBL region was of less importance as compared to 
the subtraction term). 

Finally, we comment briefly on the factorization scale 
dependence of the polarized amplitudes. We proceed in 
an identical manner to the unpolarized case in varying 
li 2 . For small £, we find a similar variation as in the 
unpolarized case for the quark singlet and larger varia- 
tions of 35 — 40% at the input scale and 1 — 30% at large 
Q 2 = 100 GeV 2 , respectively, for the gluon. At large £, 
we find variations in the real part of the quark singlet 
around 100% at the input and at large Q 2 . The imagi- 
nary part only varies between 5% at the input scale to 
20% at the evolved scale. Both real and imaginary parts 
of the gluon amplitude at large ( are very well behaved 
and vary only between 5—20% at both the input scale and 
at large Q 2 . It is encouraging that the variations due to 
factorization scale changes can be safely neglected since 
the polarized distributions are even less well known than 
the unpolarized ones. 

IV. GENERAL COMMENTS ON NLO AND 
NNLO CORRECTIONS 

In this section we make some general comments about 
the structure of NLO corrections and expected NNLO 
corrections. 

In |l^] we pointed out that the relative shape change in 
the NLO evolution of GPDs, compared with LO, is due to 
a new class of integrable divergences, In (1 — X/C,) n /(l — 
X/C,) % , n,i = 0,1,2, appearing in the region around 
£. The same type of integrable divergence also ap- 
pears in the NLO coefficient functions, but is absent at 
LO, although one has an integrable singularity of the 
1/(1 — X/Q type. This fact alone helps to explain why 
one finds, in certain regions, large changes in going from 
LO to NLO in the amplitudes. 

For DVCS observables the appearance of the gluon at 
NLO also changes things dramatically since the gluon 
contribution turns out to be of the same order as the 
quarks, at least at small £ due to an extra factor of 1/C 
in eq. (|^). So, not only is a new quantity introduced, but 
one of the same order of magnitude as our LO quantities. 
Furthermore, with the unpolarized real part of the gluon 
amplitude being negative and of comparable size at small 
C as the NLO correction to the real part of the quark am- 
plitudes (which is also negative) , observables sensitive to 
the real part, such as the azimuthal angle asymmetry 
are expected to change dramatically in NLO (see [p^[). 



Conversely, the effect should not be as dramatic for ob- 
servables sensitive to the imaginary part, despite the fact 
that the NLO gluon correction is big and negative, since 
the correction to the quark singlet at NLO is big and 
positive so the corrections will cancel to a certain extent 
(see g§). 

What would one expect in NNLO ? From experience 
obtained in calculating forward coefficient functions, and 
some evolution kernels at 0(a 2 ), one would expect to 
see only the same type of integrable divergences reap- 
pearing, maybe with different powers in logs and rational 
functions, but not a new class of divergences which could 
radically alter the behavior of the amplitudes. Also, in 
contrast to moving from LO to NLO which gives the first 
gluon contributions, no new parton species appears at 
NNLO. This leads us to speculate that the NNLO order 
corrections should be mild. 

V. CONCLUSIONS 

We have presented a detailed analysis of the quark sin- 
glet and gluon contributions to the polarized and unpo- 
larized DVCS amplitudes at NLO, using NLO-evolved 
generalised partons distributions (GPDs) built from sen- 
sible input models. We have compared throughout with 
the LO results using the same input GPDs, and have 
therefore quantized the effect of the NLO corrections. 

These results are directly relevant to measurable quan- 
tities in ep — > ep7 processes at the HERA and HERMES 
experiments, and hence may be used to constrain the 
GPDs at NLO. 

The most striking feature of our results is that for a 
given Q 2 the unpolarized amplitudes exhibit an effective 
single-power behaviour in the skewedness parameter, £, 
over a very large range, apparently indicating a universal 
behaviour. 
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A. LO and NLO coefficient functions 

The coefficient functions in eq. (^) are expanded in 
powers of a s . Up to terms of 0(a 2 ), they read 

T 1,V/A {z) =T *,LO (z) + ^) Tq ,NLO,V/A^ Q 2 f ^ > 
T 9,V/A, z) = T 9,LO {z) + ^l Tg ,NLO,V/A { ^ Q 2 j ^ ^ 

(11) 

with z = 2X/C- 1. 

The LO and NLO coefficient functions were taken from 
|l9fl and are given by: 



1 



T q ' (z) 

± — a 

T 9 > LO (z) = 0, 

Tq ,NLO,V {z) = Tq ,NLO,A {z) _ ln 1 



NLO, A 



C f 



2(1 -z) 
z 3 

(i* z z 4 

rpg.NLOy f\ _ _rpg,NLO,A 



l + z 

1-2 

2 In — h 3 x 



ln — +2 ln — 



27 

4 



N F 



i 



i 



(^) 



ln — + ln -y- 



l + z 



In 



l + z 



T g,NLO,A( \ _ -^F 
W 2 



1 



Ini^ 



. Q 2 , 1-z \ ln 2 ^ 

In + n 2 j 2— 

/x 2 2 y 2(1 + z 2 



(12) 



where the LO quark coefficient is normalized in such a 
way that, in the forward limit, after properly restoring 
the dependence on both skewedness parameters, one re- 
covers the LO DIS coefficient 8(1 — x) and the NLO gluon 
coefficient is normalized such that one recovers \C® 1 . 

In the interval [0, Q, the above coefficients are strictly 
real. However, in the interval [£, 1], they split into a real 
and imaginary parts, which can be easily deduced from 
cq. ©. 



There are four different integrals to be done, which are 
regulated by adding +ie to X (this implies analytically 
continuing the log terms in the lower half plane to pick 
up —in). This generates real and imaginary parts for the 
subtraction functions, l q s , given below, for the unpo- 
larised (V) and polarised (A) cases. One has to be care- 
ful to take the appropriate sheet of the Riemann surface 
for the logarithms in order to obtain the correct imagi- 
nary parts, i.e. to use the +ie prescription consistently. 
In order to ensure the correctness of the results below, 
we cross checked them both with MATHEMATICA and 
MAPLE. 



In 
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NLO, A 



(C,Q 2 ,m 2 ) 



i-C 



vr 2 / 1 

L 2 1 - - 
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(14) 
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(tt 2 + 9 + ln C - - liT 



i-C 



ln^_ 31n fl_C\ ln2 ^l-C 



— + 9 + lnC-ln 2 



i-C 
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1-3 



(15) 



/^c^(C,Q 2 ,m 2 ) 



3Li 2 1 



In (tt 2 + 9 + 3 ln ( | In 2 > + 

y +9 + 31nC-ln 2 0-^ 



(16) 



B. LO and NLO subtraction functions 

In this section, we present the subtraction functions 
needed in our implementation of the Cauchy principal 
value prescription of eq. (]?]), i.e. the integrals 

I q ^ A (0 = £ dX T^ V ' A - 1 

l°' v ' A (Q -- 



dX T 9 > V ' A . 

V C 



2A 



(13) 



I 9 



NLO, A 



(C,Q 2 ^ 2 ) = 
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(17) 
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